
a_trape=[-1 1];
b_trape=[0.5 0.5];

a_midpoint=[-1 0 1];
b_midpoint=[0 2 0];

a_euler=[-1 1];
b_euler=[1 0];

format rat;
C_j=Compute_C_j(a_euler,b_euler)    %compare it with Example 7.109
C_j_midpoint=Compute_C_j(a_midpoint,b_midpoint)
C_j_trapezoidal=Compute_C_j(a_trape,b_trape)
